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ABSTRACT 

We present an investigation of the polar crown prominence that erupted 
on 2012 March 12. This prominence is observed at the southeast limb by 
SDO/AIA (end-on view) and displays a quasi vertical-thread structure. Bright 
U-shape/horn-like structure is observed surrounding the upper portion of the 
prominence at 171 A before the eruption and becomes more prominent during 
the eruption. The disk view of STEREO_B shows that this long prominence is 
composed of a series of vertical threads and displays a half loop-like structure 
during the eruption. We focus on the magnetic support of the prominence ver¬ 
tical threads by studying the structure and dynamics of the prominence before 
and during the eruption using observations from SDO and STEREO_B. We also 
construct a series of magnetic field models (sheared arcade model, twisted flux 
rope model, and unstable model with hyperbolic flux tube (HFT)). Various ob¬ 
servational characteristics appear to be in favor of the twisted flux rope model. 

We hnd that the flux rope supporting the prominence enters the regime of torus 
instability at the onset of the fast rise phase, and signatures of reconnection 
(post-eruption arcade, new U-shape structure, rising blobs) appear about one 
hour later. During the eruption, AIA observes dark ribbons seen in absorption 
at 171 A corresponding to the bright ribbons shown at 304 A, which might be 
caused by the erupting hlament material falling back along the newly reconhg- 
ured magnetic fields. Brightenings at the inner edge of the erupting prominence 
arcade are also observed in all AIA EUV channels, which might be caused by the 
heating due to energy released from reconnection below the rising prominence. 

Subject headings: Sun: activity — Sun: corona — Sun: coronal mass ejections 
(CMEs) — Sun: hlaments, prominences — Sun: magnetic helds 
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1. INTRODUCTION 


The solar corona contains sheared or twisted magnetic helds overlying polarity inversion 
lines (PIL) on the photosphere. The sheared/twisted helds can be observed as hlament 
channels on the disk and as coronal cavities in limb observations; solar prominences are 
located within these regions. These structures warrant investigation because of their role 
in prominence eruptions, coronal mass ejections (CMEs), and solar hares. Understanding 
the topology and evolution of the prominence/cavity magnetic held structure as well as the 
thermodynamics of the plasma within and surrounding prominences prior to the eruption is 
key to understanding the initiation of solar eruptions. 


Solar prominences are relatively cool structures embedded in the million-degree corona. 
In Ha when viewed above the solar limb, prominences appear as bright structures against 
the dark background, but when viewed as “hlaments” on the solar disk they are darker 
than their surroundings. We will use the terms “hlament” and “prominence” interchange¬ 
ably in ge neral. For more d etailed reviews o n obs e rvations and modelin g of solar prom i¬ 


nence, see 


Hiravamal (11985 h ILabrosse et al.l (120101 ): iMackav et al.l (120101) : IParentil (120141 ): 


van Ballegooiien fc Sul (120141) . When observed on the solar disk with hi gh spatial resolu ¬ 


tion, hlaments show thin thread-like structures that continually evolve (iLin et al.l 120081) . 
Recent observations with the Solar Optical Telescope on the Hinode satellite and SDO/AIA 
have revolutionized our understanding of quiescent and intermediate prominences. When 
observed above the solar limb, such prominences always sho w many thin thread- like struc¬ 


tures. In some cases the threads ar e mainly horizonta l fe.g.. lOkamoto et al.ll2007l). in other 


cases they are mainly vertical (e.g., iBerger et al.ll2008l: ISu fc van Ballegooiienll2012L and ref¬ 


erences therein). Hedgerow prominences consist of many thin vertical threads organized 
in a vertical sheet or curtain. Upward-moving p lumes and bubble s have been observed 
in between the denser, downhowing threads (e.g., iBerger et al.l 120081) . Other prominences 


consist of isolated dark columns standing vertically above the PIL, and such prominences 
often exhibit r otational motions reminicent of “to rnados” in the Earth’s atmosphere (e.g., 
Li et ah 2012 : Su et ah 2012 : Panesar et ah 2013lj. The rotational motions have been con 


hrmed using Dopple r shift measurements (e.g., iLiggett fc ZirinI 1 19841 lOrozco Suarez et ah 


2 OI 2 I: ISu et al.l 120141) . Many quiescent prominences ha ve horn-like e xtensions that protrude 


from the top of the spine in to the cavity above (e.g., iBergerl 120121: ISchmit fc Gibsonl 12013 


Su & van Bal 

egooiien 

2012) 

nence ( 

Berger 

2012 

)■ 


egooiienll2012l) . These horns may outline a flux rope located above the promi- 


The cool prominence plasma must somehow be supported against gravity because with¬ 
out such support the plasma would fall to the chromosphere on a time scale of about 10 
minutes. It has long been assumed that prominences are threaded by horizontal magnetic 
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fields and electric currents, which provide the upward Lorentz force needed to counter gravity 


fe.g., iKippenhahn fc Schliiterl 119571: iKuperus fc Raadulll974 IPneumanl 119831) . The horizon- 
tal threads may be understood as plasma being supported at th e dips of sheared arcade 


flAntiochos et af 


van Ballegooiien 


Il994 : Aulanier 


et a 


20041: Gibson et ah 


20021) or twisted flux rope flKuperus fc Raadulll974 
20061). but the structure of the vertical threads is not 


yet fully understood. For example, for tornado-like prominences it is unclear how hori- 
zontal fields can surv i ve in the presence of rotational mo tions of the vertical structures 
(iLiggett fc Zirinlll984l) . Ivan Ballegooiien fc Cranmerl (120101) proposed that hedgerow promi¬ 
nences are located in vertical current sheets, and tha t the p lasma is supported by small-scale 
“tangled” magnetic helds within the sheet. iBergerl fl2012l) proposed that the current sheet 
is lo c ated below an elevated flux rope. A three-dimensional MHD simulation performed by 
FanI (120121) suggests that the vertical threads and horn-like structure can be explained by a 
magnetic flux rope with hyperbolic flux tube (HFT), which is the g eneralization of a 2.5D 
X-point or a 3D separator. The vortex-filament concept proposed bv lSu et al.l (120141) is sim¬ 
ilar to the “tangled field model” in the sense that the magnetic fields supporting the vertical 
threads are vertical and twisted. One problem with these ideas is that small-scale tangled 
fields are not consistent w ith the apparent smoothne ss of the magnetic fields derived from 


Hanle measurements (e.g.. lOrozco Suarez et al.ll2012l ) 


Several mechanisms have been proposed to trigger solar eruptions. The initiation mech¬ 
anisms of the explosive eruption can be divided into two categories: (1) ideal MHD instability 
or loss of equilib rium: (2) non-ideal f ast magnetic reconn ection. The hrst group includes the 
kink instabilitv( Fan fc Gibson! 2007 ). toru s instability (Kliem &: Torok 2006 


Kliemetah 


2014 ) o r loss of equilibrium (jPorbes 1990 : Lin fc ForbesI 200o|): while the flux emergenc e 
model ( Ghen fc Shibata 200ol). flux cancelatio n mode hlvan Ballegooiien fc Martens 1989 ). 
and the breakout model (lAntiochos et al.l 119991) belong to the second category. At present, 
it is often unclear which mechanism, or combination of mechanisms, is responsible for any 
particular event. 

The magnetic held plays a primary role in hlament formation, stability, and eruption 
(jPriest lll989l:lTandberg-Hanssenlll995l:lMackav et al.ll2010l) . However, the magnetic structure 
of prominences is still not fully understood, with many observations and theoretical models 
differing on the exact nature of the magnetic held. In this paper, we will study an eruptive 
polar crown prominence consisting of vertical threads and examine the possibilities of three 
types of prominence support models, namely: sheared arcade model, twisted hux rope model, 
and tangled held model. The initiation mechanism of the explosive eruption will also be 
investigated. 
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2. Observations 


2.1. Data Sets and Instrnments 


A large polar crown prominence was observed 


Atmospheric Imaging Assembly (AIA, iLemen et al. 


(SDO), as well as the STEREO_B (Behind)/EUVI f Wuelser et ahll2004j: [Howard et al.lbOOSll . 


o ern pt on March 11-12 in 2012, by the 


20121) aboard the Solar Dynamics Obse rvatory 


At the time of ernption, STEREO_B is observing the Sun ahead of the Earth, and the sep¬ 
aration angle with Earth is 118°. The photosphe ric magnetic held information is provided 
by the Helioseismic and Magnetic Imager (HMI, ISchou et ahll2012l) aboard SDO. 


The apparent slow rise of the large prominence begins around 17:00 UT on 2012 March 
11. This prominence eruption is associated with a CME, which has a median velocity of 
399 km s“\ according to the CACTus SOHO/LASCO CME catalog. The hrst appearance 
of the associated CME in the LASCO/C2 held of view is at 01:25 UT on 2012 March 12. 
Two bright ribbons and post-eruption arcade are observed during the eruption, while no 
associated hare can be clearly identihed. In the current paper, we study the structure and 
dynamics of the prominence before and during the eruption. 


2.2. Structure and Dynamics of the Prominence-Cavity System before the 

Eruption 

Figure 1 presents multi-wavelength SDO/AIA observations of the target polar crown 
prominence system at 17:01 UT on 2012 March 11 prior to the eruption. This hgure shows 
that the prominence is composed of a series of vertical threads. In 131 A, 94 A, and 304 A im¬ 
ages, the prominence displays dark absorption in the middle and bright emission in the 
surroundings. The extended bright emission surrounding the dark hlament is best seen at 
171 A, which shows that a bright U-shape/horn-like structure appears to be located above 
the dark vertical prominence threads. In 335 A, 211 A, and 193 A images, the prominence is 
seen in absorption, and the dark vertical threads are surrounded by a dark region in contrast 
to the bright emission shown at 171 A. At these three channels, we also see a cloud of bright 
emission (best seen in 211 A images) located above the vertical threads within the dark 
cavity. Some of the prominence vertical threads can be seen in emission at 1600 A. 

AIA images at 304 A (left), 171 A (middle), and 193 A (right) at 12:00 UT from March 
6 to Mar 11 in 2012 are presented in Figures 2-3. These two hgures show that the height of 


^http://sidc.oma.be/cactus/catalog.php 
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prominence (304 A) appears to increase with time, and the size of the cavity (193 A) also 
appears to become bigger as time goes on. Due to the large time range, the prominence at 
the limb is a projection of different portions of the long prominence. Therefore, the apparent 
increase of the prominence height and cavity size may be due to two reasons: the height of 
different parts of the prominence is different; or the whole prominence is rising with time. 
The end view of the prominence on March 6 shows that the vertical threads are located at the 
bottom of the dark cavity (Figure 2c). A series of nearly vertical horn-like structure appears 
to go through the entire vertical threads from bottom to top on the southern side of the 
prominence in 171 A image (Figure 2b). The side view of the prominence at 171 A on March 
7-8 shows nearly horizontal bright horn-like structure protruding the top of the vertical 
prominence threads. The cavity appears to become much bigger on March 8 (Figure 2i), 
although lots of bright loops along the line of sight are projected within the cavity. As the 
Sun rotates from March 9 to March 11, AIA observes the prominence from side view to end 
view (Figure 3). The vertical threads and horn-like structure are much wider on March 9 
(Figures 3a-3b). Figures 3b-3c shows that in 193 A image bright U-shape structure appears 
at both the top and bottom parts of the vertical threads, while at the other part bright 
U-shape structure is shown at 171 A. This suggests that the U-shape structure goes through 
nearly the entire vertical threads when placing the two images together, and the same result 
can also be obtained from the observations on March 10-11. On March 11, the structure 
containing the vertical threads appears much narrower since the threads gradually align with 
each other along the line of sight, and the horn-like structure becomes more like a narrow 
U-shape structure. The size of the cavity outlined by the big bright loop is similar to that 
on March 8, though it is still hlled with smaller bright loops projected along the line of sight. 


2.3. Structure and Dynamics of the Eruption 

SDO/AIA and STEREO_B/EUVI observations of the prominence eruption (before: top, 
during: middle, after: bottom) in three wavelengths (304 A: left, 171 A: middle, 193 A: right), 
are presented in Figures 4-5, respectively. Details of the eruption process can be found from 
the online videos 1 and 2. The end view by AIA in Figure 4 shows that at 00:30 UT on March 
12 the prominence rises much higher and becomes much wider with two dark columns on the 
two sides and fewer dark thin threads in between. Small post-eruption loops are observed 
near the limb at 193 A (Figure 4f). The prominence leaves the AIA held of view (FOV) at 
04:00 UT, the 193 A image shows much brighter and larger post-eruption loops (Figure 4i), 
and a helmet streamer structure appears in 171 A image (Figure 4h). 

The top/side view by EUVLB in Figure 5 shows that this prominence is very long. 
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and the part where the prominence ernpts hrst is near the western end, which appears to 
be higher than the other part at 17:00 UT on March 11. This eruption appears to be an 
asymmetric eruption, during which only the western leg lifts off, while the eastern part 
remains attached to the surface as shown in Figure 5d (also see video 2). After the eruption, 
series of bright post-eruption loops are observed in the 195 A images (e.g.. Figure 5i). 


The kinematics of the erupting prominence is presented in Figure 6. We trace the leading 
edge of the prominence to infer its bulk motion by hrst selecting the linear slice (black) that 
best characterizes the overall trajectory. Two additional slices (blue and green) offset by 
2° in either direction are processed identically for error estimation (Figure 6a). Emission 
along these lines at a given time is binned to 300 pixels and interpolated onto a uniform 
distance grid, yielding a spatial resolution of ~2". Light curves are drawn from 10-frame 
(2-min) averages to improve signal-to-noise and are stacked against subsequent observations 
to produce height-time images like the one shown in Figure 6b for the main slice (black). The 
height-time images are then further processed to improve contrast. Each row is multiplied 
by its height to boost signal far from the limb and t he image is t hresholded above 1.5 x its 
median value. The Canny edge detection algorithm (ICa.nnvlll986l) is then applied to extract 
the leading edge. Figure 6c shows the application of the Canny algorithm to the image in 
Figure 6b, and the pixels highlighted in red are used as the individual height measurements. 
These points are extracted automatically, but their time range must be selected manually, 
which is particularly important for the start time because it effecti vely dehnes the onse t of 
the slow-rise phase. The aforementioned procedure is also used by iReeves et al.l (120151) for 
IRIS observations of an eruptive prominence and bv lMcCauIev et ah (1201511 for a statistical 
study that includes this event. 


The height m easurements are then ht with an analytic approximation presented by 


Cheng et al.l (120131 1 for their study of an active region flux rope eruption: 


h{t) = coe^^ + ci{t- to) + C 2 (1) 

h(t) is height, t is time, and r, to, Cq, Ci, and C 2 are free parameters. This model combines 
a linear equation to treat the slow-rise phase and an exponential to treat the fast-rise. The 
onset of the fast-rise phase can be dehned as the point at which the exponential component 
of the velocity equals the linear (i.e., the total velocity equals twice the initial), which occurs 
at: 


Anset = rln(cir/co) + to 


(2) 
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Fitti ng is accomplishe d using MPFIT, a non-linear least squares curve fitting package 


for IDL flMarkwardt 


20091). Figures 6d, 6e, and Figure 6f show the £t result and its time 


derivatives. Based on this approximation, we find that the initial slow-rise velocity is 2.6 ± 
0.2 km s“^, and the maximum velocity in the AIA FOV is 110 ± 5 km The onset of 
the fast rise phase occurs at 22:57 UT ± 7 minutes at a height of 97 ± 5 Mm. At the onset 
point, the acceleration is 1.2 ± 0.1 m s“^, and the hnal acceleration is 49 ± 5 m s“^. 

Two strategies are employed to quantify uncertainty, and the errors quoted are the sum 
of both. The hrst is to identically process two adjacent slices offset by 2° on either side 
of the original. Standard deviations from these results account for ~55% of the velocity, 
acceleration, and time uncertainties. The second strategy is to perform 100 Monte Carlo 
(MC) simulations to estimate uncertainties from the £t parameters by randomly varying our 
height measuremen ts within some assumed error and rehtting Equation [H as was done by 
Cheng et al.l (1201311 . Since there are no standard errors for height measurements obtained 
by our edge detection method, the assumed height errors are chosen to yield a reduced chi- 
squared {x‘r) of 1.0 for the £t. That gives us an uncertainty of 4.1", which corresponds to 
~7 AIA pixels and ~2 pixels on the height-time image. This value is included in the onset 
height uncertainty, to which it contributes 65%, while uncertainties from the 3 separate slices 
and MC realizations account for 15% and 20%, respectively. 

In order to test the accuracy of our computed trajectory, we calculate the expected ar¬ 
rival time of the CME in SOIIO/LASCO/C2 using two different assumptions. If we assume 
a constant acceleration of 49 m s“^ after leaving the AIA FOV, the prominence would arrive 
at the LASCO/C2 FOV at around 02:14 UT on March 12th with a velocity of ~280 km 
s“^. If, instead, the prominence continues along our height-time £t, it would appear in C2 
at 02:05 UT with a velocity of ~420 km s“^. The actual arrival time of the prominence 
at 2.5 Rq is between 02:00 and 02:12 UT, and the velocity listed in the CACTus CME 
catalog is 399 ± 65 km s“^. The CACTus velocity is derived from linear hts to measure- 
ments across the entire CME from both the C2 and C 3 detectors, which cover 2.5 to 30 R© 
flRobbrecht fc Berghmansll2004l: iRobbrecht et al.ll2009fl . The prominence may then have fol¬ 
lowed our projected trajectory until nearly 2.5 Rq, after which it transitioned to a constant 
velocity or decelerated somewhat at large heights, but this is speculative given the lack of 
observations between AIA and LASCO/C2. 


2.4. Signature of Reconnection: U-shape and Blobs 

Figure 7 shows rotated views of AIA images from 19:00 UT on March 11 to 00:39 UT 
March 12, and details can be found at corresponding online video 3. The hrst and second 















rows show radial-filtered and running-difference images at 171 A, and following the same 
format AIA images at 193 A are presented at the third and fourth rows. The prominence 
gradually rises up from 19:00 UT to 23:00 UT on March 11. A post-eruption arcade (blue 
arrow) begins to appear around 23:50 UT on March 11. A newly formed U-shape structure 
(green arrow) begins to show up at the lower part of the rising vertical threads around 00:10 
UT on March 12. 


Similar to Figure 7, Figure 8 shows rotated AIA images from 00:43 UT to 01:19 UT on 
March 12. Small bright blobs (white arrow) begin to appear at the bottom of the erupting 
prominence vertical threads around 00:51 UT. The blobs are best observed at 171 A, and 
the bigger and brighter blobs can also be seen at 193 A. These bright blobs then rise up into 
the less dense prominence region located between the two dark dense vertical threads. The 
blobs appear to rise faster than the bulk motion of the erupting prominence (see video 3). 

A detailed kinematics study of the rising blobs is presented in Figure 9. Figure 9a shows 
AIA 171 A running-difference image at 01:20 UT on 2012 March 12. Distance-time plots of 
emission at 171 A, 193 A, 211 A along the white slice marked in Figure 9a are presented 
in Figures 9b, 9c, 9d, respectively. We perform linear fits (marked as red lines) to three of 
the rising blobs for each channel, and the velocity for each rising blobs is presented in the 
figure. This figure confirms that the speed of these bright rising blobs is larger than that of 
the leading edge and bulk motion of the prominence. The rising blobs that appear earlier 
(95 km s“^) are slower than those appear later (259 km s“^). The appearance height of these 
blobs is increasing with time as shown in Figures 9b-9d. In addition, bright falling features 
are also observed starting around 01:10 UT (at 171 A) on March 12 as shown in Figure 9b. 


The concave up U-shape structure, rising blobs, and post-eruption arcade are all signa¬ 
tures of reconnection. The aforementioned observations suggest that magnetic reconnection 
begins around 23:50 UT on March 11 near the bottom of the vertical threads. The increase 
of the appearance height of bright blobs may be due to the increase in height of the recon¬ 
nection point. Intermittent plasmoids or blobs are generally explained by the tearing-mode 
instability of the thin current sheet whe r e a series of magne tic islands are recurrently created 
during reconnections flFurth et al.lll963l: iDrake et al.ll2006l) . Small plasmoids or magnetic is¬ 
lands have be en found to flow along the current sheet either sunward or anti-sunward in both 


observati ons fISavage et a. 


ulations (IShen et ah 


2011 


2010 : LiuboiS : Liu et ah 2013 . and references therein) and sim- 
Karpen et ah [ESI and references therein). The speed of rising 


blobs in our eruption is at the lower end of the anti-sunward blobs’ speed in the literature 
(100-1400 km s“^). This is not surprising, since our event is a slow eruption that occurred 
at the polar crown. 
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2.5. Darkenings and Brightenings during the Eruption 


Figure 10 shows AIA observations of bright and dark ribbons during the eruption. In 
304 A (Figure 10a) image, we can see two bright ribbons located at the footpoints of the 
post-eruption arcade, which can be seen most clearly at 193 A (Figure 10c). Corresponding 
to these bright ribbons, two dark ribbons are seen in absorption at 171 A (Figure 10b). 
Distribution of intensity versus time along a slice (from south west to north east) nearly 
perpendicular to the ribbons (white line in the top left corner image) at 304 A and 171 A are 
presented in Figures lOd-lOe. Figure lOd shows that the two bright ribbons hrst appear 
around 01:30 UT, and the two ribbons then gradually move away from each other. The slope 
of the outer edge of the ribbons, i.e., the newly formed ribbons, shows that the southern 
ribbon appears to move faster. The southern dark ribbon seen in absorption at 171 A hrst 
appears around 01:10 UT, and the appearance of the full northern dark ribbon is around 
01:30 UT. At 171 A, the ribbons are only seen in absorption initially, but later on bright 
emission also shows up. Most of the northern ribbons are seen in absorption, while only the 
outer edge of the southern ribbon, i.e., the newly formed ribbons, are dark. 


The observational characteristics of dark ribbons at 171 A suggests that they might 
be caused by the erupting hlament material falling back along the newly reconhgured mag¬ 
netic helds. The co-spatial brightenings at 304 A might be partially caused by heating of 
the plasma due to the kinetic energy of falling hlament material compressing the plasma 
flGilbert et al.ll2013l) . In addition, thermal and/or non-thermal energy released from recon¬ 
nection impacting the chromosphere might also contribute to the observed brightenings, as 
suggested by the lateral appearance of brightenings at 171 A. 


Figure 11 shows brightenings of prominence material observed by AIA (left: 211 A, 
middle: 304 A) and EUVLB (right: 304 A) during eruption. The top row shows images 
before the appearance of the brightenings, and the brightenings are shown in the bottom 
row. Brightening B1 hrst appears around 01:15 UT, and brightening B2 shows up (01:19 
UT) immediately after. Both brightenings are observed in all EUV channels, and B1 appears 
to be brighter than B2, especially in the hotter channels (e.g., 94 A, 335 A). The change 
of the prominence seen from absorption to emission (B1 in Figures 11a and lid) suggests 
that the brightenings might be caused by heating. EUVLB observation shows that these 
brightenings occur at the inner edge of the erupting prominence. Therefore, the observed 
prominence brightenings are likely to be caused by heating due to energy released from 
reconnection below the rising prominence. 
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3. Modeling 


3.1. Flux Rope Insertion Method 


Ou r magnetic field models are constructed using the flux rope insertion method devel¬ 
oped bv ivan BallegooiienI (120041) . Detailed description of the methodology c an be found in the 


literature fiBobra et al 


2008 


Su et al.l I 2 OIII: ISu fc van BallegooiienI 120121) , and we describe 


the method briefly below. First, the potential held is computed from the high-resolution 
(HIRES) and global magnetic maps. Then, by appropriate modifications of the vector po¬ 
tentials a “cavity” is created above the selected path, and a thin flux bundle (representing 
the axial flux of the flux rope ($ 0 x 1 )) is inserted into the cavity. Circular loops are added 
around the flux bundle to represent the poloidal flux of the flux rope {Fpoi). The above held 
configuration is not in force-free equilibrium. So our next step is to use magneto-frictional re¬ 
laxation to evolve the held toward a f orce-free state. This method is an iterative relaxation 
method (Ivan Ballegooiien et al. 2000 ) specifically designed for use with vector potentials. 
Specifically, we solve the following equation: 


dA B 

— =r]oV X B - 772 V X B + —V ■ + V{r]dV ■ A), 


( 3 ) 


where v is the plasma velocity, rjo, 772 , and rjd are constants in space, and a = j ■ B/B'^, 
where j = V x B. The velocity is given by 


V = (fj — vif X B) X B / B‘‘ 


( 4 ) 


where / is the coefficient of magnetof riction, and vi des cribes the effects of buoyancy and 
pressure gradients in the photosphere flBobra et al.l 120081) . Magnetofriction has the effect of 
expanding the flux rope until its magnetic pressure balances the magnetic tension applied 
by the surrounding potential arcade. Significant magnetic reconnection between the inserted 
flux rope and the ambient flux may occur during the relaxation process. Therefore, the end 
points of the flux rope in the relaxed model may be different from that in the original model. 


The lower boundary condition for the HIRES region is derived from line-of-sight (LOS) 
photospheric magnetograms obtained with the SDO/HMI. Since the prominence is observed 
near the east limb, we use magnetograms that are taken several days after the prominence 
eruption on March 12. We combine four magnetograms taken on 2012 March 16-19 (each at 
18:11 UT) to construct a high-resolution map of the radial component B^ of the magnetic 
held as a function of longitude and latitude at the lower boundary of the HIRES region. We 
also use a SDO/HMI synoptic map of B^ to compute a low-resolution global potential held, 
which provides the side boundary conditions for the HIRES domain, and also allows us to 
trace held lines that pass through the side boundaries of the HIRES region. 
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The HIRES magnetic map is shown in Figure 12. Note that this region is at the polar 
crown, and the held has mixed polarity with dominantly positive polarity on the south side 
of the PIL and negative polarity on the north side. Based on this map alone, it is difficult to 
recognize exactly where the PIL is located. Therefore we used the hlament channel observed 
by STEREO_B to locate the base of the prominence/hlament on the magnetic map. The 
blue curve is the path along which the hux rope will be inserted into the model. At the 
two ends of the path (blue circles), the flux rope is anchored in the photosphere. Because 
this region is very close to the South pole, strong numerical artifacts begin to appear at 
the side near the South pole during the magneto-frictional relaxation process. To reduce the 
artifacts, we adopt a special relaxation procedure. During the hrst 70000-iteration relaxation 
we move the HIRES region close to the disk center by setting the latitude of the center of the 
HIRES box to be zero. According to our previous experience we hnd that it normally takes 
about 70000-iteration relaxation for the magnetic helds in the quiescent polar crown hlament 
system to approach a force-free state. Then we move the HIRES region back to the original 
location, then run 10000 more iterations to make the helds to interact and merge with the 
real global polar magnetic helds. After this extra 10000-iteration relaxation, the magnetic 
helds in the HIRES region become adjusted to the real global polar helds in the surrounding 
low-resolution region. If we continue to relax the magnetic helds, numerical artifacts begin to 
appear. The two criteria that we use to choose the specihc iteration numbers are: minimizing 
numerical artifacts and keeping the HIRES region being adjusted to the original surrounding 
polar magnetic helds. The parameters used during the relaxation process are presented in 
Table 1. 

Table 1: Parameters of Relaxation Method._ 


Relaxation 

Iteration 

Vo 

V2 

V4 

0-100 

1 

0 

0 

100-1000 

1 

0 

0.003 

1000-10000 

1 

0 

0.001 

10000-20000 

1 

0 

0.0003 

20000-70000 

1 

0 

0.0001 

70000-80000 

1 

0.003 

0.01 
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3.2. Models versus Observations 


We construct five models with different combinations of axial and poloidal fluxes of 
the inserted flux rope. In all cases the inserted flux ropes have sinistral orientation of the 
axial held. Table 2 shows the model parameters of the hve models. The initial inserted hux 
bundles of Models 1 and 2 are straight (untwisted) bundles, while for Models 3-5 the initial 
in serted hux bund els are twisted hux ropes with nonzero poloidal hux (Fpo^). As mentioned 


Su et al.l (1201111 . one constraint on the stability of the models is that the magnetic energy 


m 

of the held after the relaxation should be less than that of the open held. To estimate 
the energy of the open held, we change the negative polarity helds of the HIRES region to 
positive polarity and then compute a potential held. The resulting open-held model has a 
magnetic energy Eopen = 9.77 x 10^^ erg, whereas the standard potential held has energy 
of Epot = 1.15 X 10^^ erg. Therefore, the free energy of the open held is about 8.62 x 10^^ 
erg. This requires that the free energy of the hux-rope models be less than 8.62 x 10^^ erg. 
Note that these energies refer to the HIRES part of the computational domain only (not the 
whole Sun). All of the hve models presented in this paper meet this criterion. 


Figure 13 presents the current distribution in a vertical plane indicated by the yellow 
line (nearly perpendicular to the PIT) in Figure 12 for Models 1-4 after relaxing for 80000 
iterations. The strong current is concentrated in the white region. The black and white 
vectors refer to the magnetic vectors. This hgure shows that Model 1 is close to a normal 
sheared arcade, while Model 2 appears to be a hux rope conhguration with an X-point/HFT 
and sheared arcade located below. Model 3 and Model 4 display an elevated twisted hux 
rope conhguration. 

Figure 14 shows comparisons of Models 1-4 with SDO/AIA and STEREO_B observa¬ 
tions. The colored curves in Figure 14 show held lines selected by clicking diherent points 
in the two dimensional plots shown in Figure 13 from the four models. The blue features 


Table 2. Model Parameters for the Polar Crown Filament System. 


Model 

No. 

^axi 

(10^0 Mx) 

^pol 

(10^*^ Mx cm“^) 

^ free 

ergs) 

Afreet Epoten 
(%) 

1 

2 

0 

1.98 

17 

2 

3 

0 

2.98 

26 

3 

2 

2 

8.14 

71 

4 

3 

3 

16.2 

141 

5 

3 

1 

5.76 

50 


^The potential field energy is 1.15^^ erg 
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indicate dips in the field lines (depth of color increases with height). The red and green con¬ 
tours refer to the HMI photospheric flux distribution in the HIRES region. The background 
images in the row 1, row 2, and row 3 are taken at 171 A, 193 A by SDO/AIA, and at 
171 A by STEREO_B, respectively. We can see that the dips of the field lines in Model 1 is 
lower than the prominence observed by both SDO/AIA and STEREO_B/EUVI. Moreover, 
the cavity in the model is much smaller than observed, and the model cannot reproduce the 
observed U-shape/horn-like structure. Therefore, we think that the sheared arcade topology 
does not match the observations. For Model 2, the height of the dips of field lines matches 
the AIA observations well, but is much lower than the height of prominence seen in the 
STEREO_B observations, especially the part wher e it erupts first. This model exhibits an 
X-point topology as suggested by the simulation in iFanI (120121) . and the location and height 
of the U-shape structure in the model match the SDO/AIA observations well, but not the 
STEREO_B observation. Moreover, the size of the “U” in the model is much smaller than 
that of the observation. The height of the dips of field lines in Models 3 and 4 correspond 
very well with the height of the observed prominence by AIA and EUVI, though a small 
offset towards the south in the models is identified when compared with AIA observations. 
There appears to be different sizes of cavities in the observations, since the prominence is 
very long. The size of flux rope in Model 3 appears to match the foreground cavity, while 
for Model 4 it appears to match the background larger cavity better. 


Figure 15 shows comparisons of the “U” structure in Model 4 (right column) and Model 
5 (left column) with that in observations. The color curves in this figure are selected field 
lines that going through a vertical plane near the location of “U” structure where the promi¬ 
nence firstly erupts (white arrows in Figures 15a-b). The bright “U” structure in both AIA 
and EUVI observations are indicated with the solid white arrows. Corresponding U-shape 
structure in the models are marked with dashed white arrow (light blue and pink lines). This 
figure shows that the location, size, and height of the U-shape structure in Model 4 match 
the observations well. While for Model 5, the “U” structure is much lower than that in the 
observations. 


4. Discussions 


4.1. Onset of Fast Rise: Torus Instability 


In Section 2.3, we found that the prominence eruption beg ins with slow rise then 


evolve to fast rise , which is often seen in solar filament eruptions fjSterling fc Moorel 12005 


Cheng et al.l 1201311 . The important question is what causes the transition from slow rise to 
fast rise? In other words, what causes the onset of explosive fast rise phase? In this section. 
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we investigate two possibilities: i.e., torus instability and magnetic reconnection. 


The threshold of torus instability is given in terms of the decay index of the external 
poloidal held, at the position of the current channel, n = —dlnB/dlnz > n^r- The canoni¬ 


cal values of the critical decay index a re 1.5 for a toroidal current channel ( Bateman 19781) 
and 1.0 for a straight current channel ( van Tend fc KuperusI IT^ . The critical index {ucr) 
for to rus instability ranges from 1.0 to 2.0 in theoretical calculations and numerical simula¬ 
tions ( Kliem fc Torok 2006 : Fan fc Gibson 2007 : Aulanier et ah 2010; Demoulin fc Aulanier 
2010l) . On the other hand, from the measurement of the height of a s e t of q uiescent promi¬ 
nences, combined with potential held extrapolations, iFilippov fc Peru (120011) found Ucr ~ 1- 


Figure 16 shows distributions of the magnitude of horizontal components of the potential 
held (top) and the decay index (bottom) of the horizontal components with height over the 
main PIL at two locations marked with a yellow line (middle) and a star sign (right) in 
Figure 12. The dashed and solid vertical lines refer to the height of the outer edge of the 
prominence and the center of the dark cavity at the onset of fast rise phase, respectively. 
Due to the variety of the supporting magnetic conhguration, we take the hrst (used by 
Filippov fc Deni (120011 )1 and second heights as the lower and upper limits of the apex of the 
flux rope magnetic axis, respectively. We hnd that at the onset of fast rise phase, the decay 
index n is 1 ± 0.2, which is very close to the critical value of a straight current channel for 
torus instability. STEREO observations suggest that this prominence begins to erupt near 
the west end rather than the middle. Figures 16 shows that the lower limit of the decay 
index near the west end is slightly larger than that in the middle at the onset of fast rise. We 
should note that the potential held model is based on the magnetograms taken several days 
later, which is only a close approximation. The magnetic helds at the time of the eruption 
may be different due to disturbance of a prominence eruption nearby several hour earlier. 
Therefore, we cannot exclude the possibility that the overlying held actually falls oh more 
quickly with height as the the reason why the western portion of the hlament erupts hrstly. 


4.2. Magnetic Structure Supporting the Vertical Threads 


The aforementioned comparisons between models and observations suggest that the 
twisted hux rope model (Model 4) best matches the observations, while t he sheared arcade 
model (Model 1) shows the worst match in comparison to observations. IXia et al.l (120141) 
presents a simulation of the in situ condensation process of solar prominences. In this 
simulation the vertical prominence resides in the horizontal helds of concave upward held 
regions of a hux rope. The synthetic SDO/AIA view at 304 A and 211 A in Figure 4 by 


Xia et al.l (120141) match our observations (Figure 1) very well, but not for the other two 
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channels, i.e., 171 A and 193 A. The synthetic SDO/AIA views show a horn and cavity 
strnctnre located above the vertical prominence threads at 193 A and 211 A. While in 
observations, the “U” strnctnre is shown at 171 A. This prominence is similar to the tornado- 
l ike prorninences which are oft en observed to show rotational motion along the vertical axis 
(iLi et al.ll2012l:ISu et al.ll2012[l . It is nnclear how horizontal helds of the flnx rope can snrvive 
in the presence of rotational motions of the vertical strnctnre. 


Model 2 is closest to the “tangled held model” proposed bv ivan Ballegoohen fc Cranmer 


(l2010l) . in which the vertical prominence threads are snpported by the tangled helds in a 
vertical cnrrent sheet below a twisted hnx rope, and the “U” s trnctn re corresponds to the 


central vertical cnrrent layer in the MHD simnlation by iFanI fl2012i) . Dne to limitations 
in the magneto-frictional method, we are not able to reprodnce the vertical cnrrent sheet 
as snggested in the “tangled held model”, instead we obtain a conhgnration with newly 
reconnected arcade located below an X-point. Moreover, the magnetic conhgnration can be 
strongly distorted by the weight of prominence becanse the magnetic held at polar crown is 
very weak. Therefore, althongh Model 2 does not match the observations well, we cannot 
rnle ont the “tangled held model” for the target prominence. Since nnlike the twisted hnx 
rope model, the “tangled held model” can explain the rotational motions of tornado-like 
vertical strnctnre. 


5. Summaries and Conclusions 

In this work, we stndy the magnetic strnctnre and dynamics of a tornado-like ernptive 
polar crown prominence nsing both observations by SDO and STEREO_B and magnetic held 
modeling. Onr main hndings are snmmarized below: 

(1) STEREO_B observes that the prominence is a very long strnctnre consisting of 
series of vertical threads. AIA observes the prominence at the sonth-east limb. Prior to 
the ernption, the prominence gradnally becomes mnch higher and the cavity becomes mnch 
bigger as the Sun rotates from 2012 March 6 to March 11. A Horn-like/U-shape structure 
appears to go through nearly all the vertical threads from top (171 A) to bottom (193 A). 

(2) The slow rise of the prominence begins around 17:00 UT on March 11, it then evolves 
to fast rise around 22:57 (± 7) UT on March 11, when the height of prominence’s leading 
edge is 97 ± 5 Mm. In the AIA held of view, the maximum velocity is 110 ± 5 km and 
the hnal acceleration is 49 ± 5 m s“^. Comparing with SOHO/LASCO CME observations, 
we hnd that the prominence may have followed our height-time ht (i.e., with increasing 
acceleration) after leaving AIA held of view until nearly 2.5 R©. 
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(3) A post-eruption arcade begins to appear around 23:50 UT on March 11. A newly 
formed U-shape structure begins to show up at the lower part of the rising prominence 
threads around 00:10 UT on March 12. Around 00:51 UT, we start to see small bright blobs 
at the bottom of the rising vertical threads. The speed of the rising blobs is faster than the 
leading edge and bnlk motion of the ernpting prominence. The blobs that appear earlier (95 
km s“^) are slower than those appear later (259 km s“^). These blobs can be explained by 
the tearing mode instability of a thin current sheet where a series of magnetic islands are 
recurrently created dnring magnetic reconnection. 


(4) During the eruption, AIA observes dark ribbons seen in absorption at 171 A cor¬ 
responding to the bright ribbons shown at 304 A. The observational characteristics of dark 
ribbons at 171 A suggests that they might be caused by the erupting hlament material 
falling back along the newly reconhgured magnetic helds. Dark ribbons are also reported by 


Xiao et al.l (l2015l ). who interpret it as a void region with smaller magnetic held strength and 


lower plasma density caused by magnetic held dehection dnring magnetic reconnection. The 
co-spatial brightenings at 304 A might be partially caused by heating of the plasma dne t o 
the kinetic energy of falling hlament material compressing the plasma (iGilbert et al.ll2013|) . 


(5) Brightenings at the inner edge of the ernpting prominence arcade are observed in all 
AIA EUV channels during the ernption. These brightenings might be caused by the heating 
dne to energy released from reconnection below the rising prominence. 


(6) Using the hnx rope insertion method, we constrnct a series of magnetic held mod¬ 
els (sheared arcade model, twisted hux rope model, and unstable model with HFT), then 
compare with both SDO and STEREO_B observations. Various observational characteristics 
appear to be in favor of the twisted hux rope model. However, the “tangled held model” can¬ 
not be ruled out, since it can explain the rotational motions of tornado-like vertical strnctnre 
which the twisted hux rope model cannot. 


(7) We hnd that the hnx rope snpporting the prominence enters the regime of torus 
instability at the onset of the fast rise phase, and signatnre of reconnection (post-eruption 
arcade, new U-shape strnctnre, rising blobs) appears about one honr later. This resnlt 
snggests that the transition from slow rise to fast rise phase of this prominence eruption is 
likely to be caused by the torus instability. 
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Fig. 1.— Multi-wavelength images taken by SDO/AIA at 17:00 UT on 2012 Mar 11 before 
the eruption. A color version of the hgure is also available in the electronic edition of the 
Astrophysical Journal. 
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Fig. 2.— SDO/AIA observations of the eruptive prominence before the eruption at 12:00 UT 
from Mar 6 to Mar 8 in 2012. A color version of the hgure is also available in the electronic 
edition of the Astrophysical Journal. 
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Fig. 3.— SDO/AIA observations of the eruptive prominence before the eruption from Mar 
9 to Mar 11 in 2012. A color version of the hgure is also available in the electronic edition 
of the Astrophysical Journal. 
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Fig. 4.— SDO/AIA observations of the prominence eruption on 2012 March 12. A color 
version of the hgure and a video of this event (video 1) are also available in the electronic 
edition of the Astrophysical Journal. 
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Fig. 5.— STEREO_B/EUVI observations of the prominence ernption on 2012 March 12. A 
color version of the hgnre and a video (video 2) of this event observed by STEREO_B/EUVI 
are also available in the electronic edition of the Astrophysical Journal. 











26 


AIA 304 A at 2012/03/12 01:03 UT 



Start time = 2012/03/11 19:01:02 Start time = 2012/03/11 19:01:02 


Fig. 6.— (a) Trajectories used for tracking. The black slice is selected manually, and the 
others are offset by 2° in either direction, (b) Height-time image for the black slice in Panel 
A. The dashed line indicates the time shown in Panel A, and asterisk indicates the fast-rise 
onset point, (c) Output of the Canny edge detection algorithm applied to Panel B. The red 
pixels are used as individual height measurements, (d) Fit to the red pixels in Panel C. The 
dotted lines indicate the onset of the fast-rise phase, and the colors correspond to results 
from the different slices in Panel A. (e) Velocity and (f) acceleration prohles for the height 
prohle in Panel D. A color version of this hgure is also available in the electronic edition of 
the Astrophysical Journal. 
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Fig. 7.— Newly reconnected U-shaped loops (green arrows) and post-ernption arcade (bine 
arrows). The hrst and third rows show AIA images at 171 A and 193 A at different time 
dnring the eruption. The images are first summed then radial filter technique is applied to 
enhance the emission above the limb, and the cadence is 36 seconds. The second and fourth 
rows show AIA running-difference images at 171 A and 193 A at different times during 
the eruption. The images within 36 seconds are firstly summed then running difference is 
applied. A color version of this figure and a video (video 3) are also available in the electronic 
edition of the Astrophysical Journal. 
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Fig. 8.— Observations of the rising blobs. The layout of this hgure is the same as that in 
Figure 7. The time of the images in this hgure is later than those in Figure 7. The bright 
blobs are marked with white arrows. A color version of the hgure and a video (video 3) are 
also available in the electronic edition of the Astrophysical Journal. 
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Fig. 9.— Height-time analysis of the rising blobs from 00:00 UT to 02:00 UT on 2012 March 
12. AIA image at 171 A at 01:20 UT on 2012 Mar 12. (b)-(d) refer to position-time plots of 
emission along the white stripe shown in (a) at 171 A, 193 A, 211 A. Velocities are derived 
from linear hts to height-time plots of several blobs marked with red lines. 
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Fig. 10.— AIA observations of bright (304 A, white arrow) and dark (171 A, black arrow) 
ribbons during the eruption, (a)-(c) AIA images at 01:40-01:41 UT on 2012 March 12 at 
304 A, 171 A, and 193 A. The contours of bright ribbons at 304 A is displayed on (c). (d) 
Stack (from 00:00UT and last 3.9 hours ) plot of the emission at 304 A along the white slice 
shown in the top left Conner of the image, (e) The same as (d) but for 171 A . The black 
arrows refer to the bright and dark ribbons. 
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Fig. 11.— Prominence brightening observed by SDO and STEREO. The hrst and second 
columns show AIA images at 211 A and 304 A, and STEREO-B images at 304 A are shown in 
the third column. The top and bottom rows show images at different time. The prominence 
brightenings are marked with letter “Bl” and “B2”. A color version of the figure is also 
available in the electronic edition of the Astrophysical Journal. 








32 



-80 -60 -40 -20 0 20 

LONGITUDE (deg) 


Fig. 12.— Longitude-latitude map of the radial component of magnetic field in the pho¬ 
tosphere in the HIRES region of the model. The longitude zero-point corresponds to the 
central meridian on 2012 March 16 at 18:11 UT. The blue curve shows the path along which 
the flux rope is inserted into the model. The purpose of the yellow line and star signs will 
be described in the caption of later hgures. 
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Fig. 13.— Slices of current distributions along the yellow line in Figure 12 for four different 
magnetic field models after 80000-iteration relaxations. The initial poloidal and axial flux 
of the inserted flux is shown in each panel. 
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Fig. 14.— Comparison of the four magnetic field models with observations. The background 
images at the hrst, second, and third rows are observed at 171 A and 193 A by AIA, and 
at 171 A by STEREO-B, respectively. The color curves and blue features refer to selected 
magnetic held lines and held-line dips from four different models. The red and green contours 
show the observed HMI photospheric hux distribution as shown in Figure 12. 
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Fig. 15.— Selected field lines of flux ropes from two different models (left and right) com¬ 
pared with observations. The background images on the top and bottom panels are at 
171 A taken by AIA and STEREO-B before the eruption. The red and green contours show 
the observed HMI photospheric flux distribution as shown in Figure 12. The U-shape struc¬ 
ture in observations are marked with solid white arrow, the curved feature of the model held 
lines are marked with dashed white arrow. 
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Fig. 16.— The magnitude of total magnetic field and decay index versus height plots from 
the potential field model. The left, middle, and right plots are taken at different places which 
are marked as star signs in Figure 12. 









